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abstract 

This investigation addresses the problem of 
hydrodynamic stability of a two-component 
droplet undergoing spherically-sym metrical 
gasification. The droplet components are 
assumed to have characteristic liquid species 
diffusion times that are large relative to 
characteristic droplet surface regression 

times. The problem is formulated as a linear 
stability analysis, with a goal of predicting 
when spherically-symmetric droplet 
gasification can be expected to be 
hydrodynamically unstable from surface- 
tension gradients acting along the surface of a 
droplet which result from perturbations. It is 
found that for the conditions assumed in this 
paper (quasisteady gas phase, no initial 
droplet temperature gradients, diffusion- 
dominated gasification), surface tension 
gradients do not play a role in the stability 
characteristics. In addition, all perturbations 
arc predicted to decay such that droplets were 
hvarodynamically stable. Conditions arc 
identified, however, that deserve more 
analysis as they may lead to hydrodynamic 
instabilities driven by capillary effects. 

INTRODUCTION 

Influences of volatility differences and liquid- 
phase transport on the vaporization and 
combustion behaviors of multicomponent 
droplets have been the subject of numerous 
investigations because of their potential 
importance to the efficient use of fuel sprays 
as well as their fundamental interest. Many ol 
these investigations have focused upon 


sphericalljfa'vmmetric low Reynolds-number 
evaporation and combustion of 
multicomponent droplets. These 
investigations have demonstrated that liquid- 
phase species transport may be important 
when surface regression is fast compared to 
species transport, and when significant 
volatility differences exist between liquid 
components. Landis and Mills [1], for 
example, report results of numerical analyses 
concerning spherically-symmetric 
evaporation of miscible binarv droplets, 
demonstrating the effects of volatility 
differences and slow liquid-phase species 
transport. Law et al. have presented 
numerical and theoretical analyses for 
spherically-symmetric binary fuel evaporation 
with or without combustion [2*4]. 
Experimental work in this area has been 
active [see 5-12 for example]. Law el al. 
[12 13] have also reported results regarding 
evaporation of methanol or other liquid 
droplets in humid air; a focus of these studies 
was the effects of water condensation on 
droplet evaporation. Reduced-gravity 
experiments on methanol droplet combusuon 
where combustion-generated water was 
reasoned to have condensed onto the liquid 
droplet have been reported by Choi et al. 
[14]. 

Some of ihese investigations have 
demonstrated, for binary droplets with 
sufficiently different liquid component 
volatilities, that an initial transient period 
exists during which the more volatile 
component is preferentially evaporated. 
During this time, as the droplet volume 
decreases, the fraction of the less volatile 
component in the droplet is increased relative 
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to the initial state. Surface regression, when 
sufficiently rapid compared to liquid species 
transport, causes a "concentration boundary 
layer to develop near the droplet surface, 
where spedies profiles exhibit sharp spatial 
changes; the less volatile component is 
"swept up" In a thin layer because It cannot 
diffuse away from the droplet surface rapidly 
enough. It has been observed that under 
appropriate conditions (e.g., with species 
transport sufficiently slow compared to 
surface regression rates), a state may be 
attained where change**. in droplet 
compositions occur slowly with respect to 
changes id droplet size [1,2]. This stale, 
which occqrs after the initial transients have 
decayed, typically exhibits a concentration 
boundary-layer structure.* For sufficiently 
small diffusion rates, this boundary-layer 
structure Imay persist for most of the 
remaining droplet lifetime. 

Based on this observation, Law and Law 
[11] presented an asymptotic analysis 
concerning spherically-symmetric 
vaporization and combustion behaviors pf 
multicomponent droplets for the case where 
liquid-phase species diffusion is slow 
compared to droplet surface regression rates. 
In their analysis, Law and Law focused upon 
liquid-phase behavior for the case where the 
boundary-layer structure changes slowly, 
following decay of Initial transients. A 
conclusion of this work is that once the 
boundary-layer composition profiles have 
been established in a droplet, the mass-flux 
■ fraction of ahy species oft the droplet surface 
is equal to the initial llquid-pnase mass 
fraction of that species prior to vaporization. 
The analysis of Law and Law [11] was 
extended by Shaw [15], who demonstrated 
that following decay of initial transients, the 
mass-flux fraction of any species off the 
droplet surface is instead approximately equal 
to tne volume-averaged mass fraction of that 
species in the droplet interior. 

If temperature or species concentration 
gradients exist at the surface of a liquid 
droplet (e.g., because of small perturbations 
in the environment or within the droplet), it is 


* Here, "boundary layer" refers to the species 

concentration boundary layer near the droplet surface. 
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possible for surface-tension gradients to 
exist. If sufficiently intense, these surface- 
tension gradients can cause significant 
motions of the liquid-gas interface, 
influencing transport of species and energy in 
the gas and liquid phases. It is capillary 
effects such as these that this paper is 
concerned with. Theoretical studies will be 
performed to determine the influence of 
capillary convection on droplet hydrodynamic 
stability, with a goal of determining how 
capillary convection may compromise the 
attainment of spherical symmetry. The goal 
of this research is to perform a linear stability 
analysis to determine whether the spherically- 
symmetrical state is stable with respect to 
small perturbations in the flow variables. 
Viscous effects as well as the effects of 
surface tension gradients will be allowed. 
The analysis will consider only the portion of 
the droplet history where the spnerically- 
symmetrical species profiles may be 
considered quasi-steady. Earlier portions of 
the droplet history associated with transient 
liquid healing and the initial buildup of the 
species concentration profiles will not be 
addressed. 

It is noted that there have been many 
theoretical and experimental studies dealing 
with the motions of droplets and bubbles 
resulting from capillary effects. A literature 
review has revealed, however, that most 
previous studies of droplet gasification have 
generally neglected effects of surface-tension 
gradients, though a few studies have 
appeared. Higuera and LifltSn [16] have 
presented a linear stability analysis of an 
unsupported and stationary droplet 
vaponzing In a hot stagnant atmosphere. Chai 
et al. [17] have also presented results 
pertaining to stability of evaporating droplets 
vaporizing in a stagnant environment. 
Lozinski and Matafon [18] presented 
analytical studies of spinning droplets with 
thermocapillary effects. Niazmand et al. [19] 
have recently presented numerical 
calculations for vaporizing droplets that are 
moving relative to their environment. The 
calculations of Niazmand et al. showed that 
capillary flows can significantly affect droplet 
internal velocity and temperature fields of 
vaporizing octane and methanol droplets, 
even when droplet initial Reynolds numbers 
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(relative to the gas phase) are as large as 50. 
It is noted that following the initial 
publications of Niazmand et al. [19], Shih 
and Megaridis [20] presented computational 
results on the effects of thermocaplllary 
stresses and flows for nearly the same 
conditions as those considered by Niazmand 
et al. None of these previous studies have 
considered droplets with more than liquid 
species, which is the primary subject of this 
paper. 


ESTIMATES OF THE IMPORTANCE OF 
CAPILLARY EFFECTS 

Estimates of the importance of capillary 
effects may be made by considering a droplet 
of radius R. A temperature or species 
concentration gradient along the liquid 
surface will produce a surface-tension 
gradient that will be balanced by surface 
shear stresses in the gas and liquid phases. 
For these estimates, gas-phase shear stresses 
will be assOmed negligible relative to liquid 
shear stresses such that surface-tension 
gradients are primarily balanced by liquid 
shear. Neglect of gas-phase shear in surface 
force balances is justified if the gas-phase 
viscosity Is small relative to the liquid 
viscosity, a condition that is generally 
satisfied for subcritlcal conditions, and if 
liquid velocities relative to the ambient gas are 
small, i.e., there is not a thin momentum 
boundary layer in the gas phase. Liquid- 
' phase shear will be assumed to be 
characterized by Ulu/R, where |1 l is the 
liquid viscosity, u a velocity characteristic of 
the flow induced in the droplet and u/R 
represents a characteristic velocity gradient al 
the droplet surface. Use of R as a 
characteristic distance should be appropriate 
since for hydrocarbon and alcohol fuels of 
interest, the ratio of the evaporation constant 
k (based on the droplet diameter, that is, k = - 
4d(R 2 )/dt, where t is time) to the liquid 
kinematic viscosity v L is of order unity or 

less (i.e., k/VL^ 1)* Physically, if k/VLi 1, 
momentum changes at the droplet surface 
should have sufficient time to influence the 
droplet interior before the droplet evaporates. 
A balance between liquid snear stress and 


surface tension gradients gives u - Ao/[1l m 

r/tc where Aa is a characteristic difference in 
surface tension from one side of the droplet 
to the other, and t^ is a characteristic time for 
. a fluid particle to be convected halfway 
,s. around the droplet periphery (surface-tension 

gradients are characterized by Aa/R, and gas- 
phase shear is neglected assuming that the 
liquid-phase viscosity is much larger than the 
gas-phase viscosity). If we define tj, = R 2 /k 
as a characteristic gasification ^rijne of a 
droplet, we may expect that capillary effects 

are negligible if tjj/t c = RAa/(k[tL) « L 

Conversely, if RAa/(kp,L) ^ 1, capillary 
effects may cause significant convective 
flows over die lifetime of a droplet 

The magnitude of Aa depends upon gradients 
of temperatures and species concentrations 
along the droplet surface. Let us first 
consider surface temperature gradients alone. 

For this case we may write Aa - |9a/3T I 
AT, where T is temperature, a surface 

tension, and AT is the temperature difference 
from one side of the droplet to the other. If 
this relation is used in the expression for tb/tc, 

we may solve for AT ~ ktiL/(R l3a/9T I) as 

being required for tb/tc m L For a 1 mm 
diameter hydrocarbon or alcohol droplet 
undergoing combustion, estimates suggest 

that AT «- 0.01 K is required to make tb/tc - 
1. These estimates are based upon the 
characteristic values k = 10‘ 6 m 2 /s, Pl - 4 x 
l<Hkg/(m s), and |9 ct/ 3T I = 10* 4 N/(m 6 K), 
which are appropriate for near-atmospheric 
pressure conditions and droplet gasification 
in a hot environment or combustion. Since 

the ratio tb/t c scales linearly with AT, 

increases in AT produce increases in tb/tc. It 
thus appears that even very small temperature 
differences may induce significant convective 
flows. 

We may also consider the effects of species 
concentration gradients along the droplet 
surface, neglecting temperature gradients. 
For simplicity, attention will be restricted to 
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droplets composed of two miscible liquid 
components. For this case, we may write Act 

«* |9a/9xj |Axi where xi is the mole fraction 
of component "one" at the liquid surface, and 
Axj is the difference in xi from one side of 
the droplet to the other. Use of this result in 
the expression for tb/t c yields Axi ~ 
kHi/(R idc/dxi I) for tb/tc ~ 1 • 

Consider a methanol droplet of about 1 mm 
diameter. It is known that methanol (or other 
alcohol) droplets may absorb water from the 

gas phase. The characteristic values l9o/9xj I 
= 0.03 N/m, k « 1(H m 2 /s, and |1 l = 4 x 
1(H kg/(m s) are appropriate for making 
order-of-mbgnltude estimates for methanol- 
water mixture droplets undergoing 
combustion at near atmospheric pressure. 
Use of these values provides the result that 

Axi - 10* 5 for tb/tc •* 1. Hence, very small 
variations in water concentration along the 
surface of a methanol droplet may induce 
significant convective flows. Since the ratio 

tb/tc scales linearly with Axi, larger values of 

Axj produce increases in tb/tc- Significant 
capillary flows will negate the possibility of 
attaining liquid-phase spherical symmetry by 
enhancing mixing in an asymmetric and 
unknown Manner. Estimates similar to those 
presented above for methanol-water mixtures 
have also been made for miscible binary 
hydrocarbon mixture droplets. These 
estimates suggest that for droplet components 
that vary significantly in volatility, e.g., 
heptane-hexadecane, small differences in 

surface mole fractions of one component (Ax 

- HT 4 ) may cause tb/tc - 1- Such droplets 
(with significant volatility differentials) have 
been the subject of a number of reduced- 
gravity experimental efforts which are 
striving for spherical symmetry. 

Based upon the above estimates, it appears 
that small perturbations in temperature or 
species profiles along the surface of a binary 
miscible droplet may induce significant 
convective flows. If sufficiently intense, 
these convective flows will seriously 


compromise spherical symmetry, making 
inteipretation of reduced-gravity experimental 
results (which are presumed to be spherically 
symmetrical) difficult. For these reasons, 
theoretical studies of the influences of 
capillary phenomena on droplet gasification 
and combustion are presented below. These 
studies will allow the influences of capillary 
phenomena to be delineated with greater 
precision, allowing for better Interpretation of 
reduced-gravity experimental results, and 
also for increased understanding of droplet 
versification phenomena in general. For 
simplicity, the analyses are restricted to 
evaporation in a hot environment; effects of 
gas-phase combustion will be considered at a 
fater time. However, it is noted that the 
results presented here should be valid for 
combustion of droplets when the flame lies a 
significant distance away from the droplet 
surface, for example when droplets arc 
burned in environments with low oxidizer 
concentrations. 

In the following, the basic unperturbed 
solution for quasisteady droplet evaporation 
in a hot environment will first be presented. 
This solution, which is presented for the 
convenience of the reader, is a summary of 
an analysis presented by Shaw [15]. After 
presentation of this solution, the viscous- 
flow governing equations will be perturbed 
and linearized. These equations will be 
solved and stability criteria will be developed. 


TUB UNPERTURBED SOLUTION 


Liquid-Phase Analysis 

With the assumption of a constant liquid 
species diffusivily Dl and also assuming that 
liquid densities do not vary with lime, liquid- 
phase species conservation may be taken to 
be described by the partial differential 
equation and boundary conditions (where yi 
represents the mass fraction of species i) 


r 2 dr l dr/ 


( 1 ) 


4 


FROM U C D 11 P M E E M 6 R G 


01 . 03 . 1995 


13 : 18 


P. 6 



and the initial condition yi(r,0)=yio- Here, yio 
is taken for simplicity to be constant. The 

variable 8t Is the mass-flux fraction of species 

1 from the droplet surface i.e., 8j = 
(dmi/dt)/(dm/dt), where m is the droplet 
mass, mi the mass of component i in the 
droplet, and t the time. The droplet surface is 
regressing, so th“ location of the liquid-phase 
boundary is time dependent. In the 
following, the droplet surface is fixed at unity 

with the transformation x=r/R. Then with e=- 
DL/[R(dR/dt)], liquid phase species 
conservation is described by the initial value 
problem . 



yi(x,0) = yio . 

The physical time t and the nondimenslonal 

time t are related by t= f dx/[dln(Ro/R)/dl]. 

Jo 

When d-square law surface regression holds, 
i.e., where the droplet surface area decreases 

linearly with time, e=8DL/k, where k=- 
8R(dR/dt) is the evaporation constant. For 
this case it then follows that t/tb=l-e- 2x , 

where tb=4Ro/k is the droplet evaporation 
lifetime and Bio is the initial droplet radius. 

The situation of interest here is when liquid- 
phase species transport is slow relative to 
droplet surface regression. Consider the 
second component to have a volatility that is 
low relative to the first component. The more 
volatile component will be initially 
preferentially vaporized, leaving the less 
volatile component behind in the droplet. The 
droplet will become richer in the less volatile 
component relative to the initial state, with 
most of the lower-volatility component left 


behind being collected within a boundary 
layer at the droplet surface. As surface 
regression proceeds, the surface mass 
fraction of the less volatile component builds, 
and a state will eventually be attained where 
vaporization of the less volatile component 
becomes significant. Since both species will 
then vaporize vigorously, and liquid-phase 
diffusion is slow, further changes in droplet 
species profiles will occur slowly. Only that 
part of the combustion history where both 
components are vaporizing vigorously with 
slowly-changing species profiles will be 
analyzed here. Analyses of transient 
boundary layer buildup have been given 
elsewhere (21]. 

Inspection of Eq. (2) suggests that once 
boundary layer growth has slowed 
considerably and surface conditions change 
slowly, further changes in species profiles 
over die droplet radius will occur over times 

At which may be of order (1/e). The time 

period At »* 0(l/e) may comprise the major 
portion of the remaining droplet lifetime in 
physical time, since e«1 is expected, as 
argued below. Because characteristic changes 
in y are of order unity or less, the magnitude 

of the partial derivative dyjdx for these times 

may be estimated as (8yi/3x) ^ 0(e). Hence, 
it may be assumed that the boundary-layer 
structure which has evolved once vigorous 
vaporization of both components begins will 

be adequate for evaluating 8j, These 
observations are qualitatively supported by 
available numerical results [2], where the 
boundary-layer structures were observed to 
pereist for most of the droplet lifetime, until 
near the end, when the more volatile 
component was vaporized much faster than 
the less volatile component 

Analysis of the boundary layer structure 
proceeds by considering Eq. (2). Because Dl 
is generally small compared to characteristic 
magnitudes of R(dR/dt), it may be considered 

that £ « 1. The "stretched" coordinate z=(l- 

x)/e may then be inserted into Eq. (2). By 
using techniques of matched asymptotic 
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expansions, it may be shown [15] that the 
mass-flux fraction may be expressed as 

8i = y, 0 + 3e(yi. - yio) + 0(fi 2 ) (3) 

.while the asymptotic boundary-layer structure 
is provided by 

yi « yio + 

(y 4 . - yio)exp(-z)[l + e(3z + z 2 /2) + ...]. (4) 

To complete the formulation of the model, the 
gas phase and the gas/liquld coupling must be 
analyzed, as presented below. 


Gas Phase Analysis 

For analysis of the gas phase, the 
simplifications of unity Lewis number, 
quasisteadiness i.e., neglect of the transient 
terms in the governing equations), negligible 
body, forces, negligible radiant heat flux, 
equal and constant specific heats and binary 
diffusion coefficients, and that diffusion 
obeys Fick's law will be employed. In the 
following, the subscript + denotes that 
conditions are evaluated on the gas side of the 
gas-liquid interface. 

Assuming spherical symmetry, the energy 
and specieS conservation equations for the 
gas phase can respectively be written as 



where p is density, Cp is specific heat, X is 
thermal conductivity (assumed constant), T is 
temperature, r is the radial distance outward 
from the droplet center, D g is the gas-phase 
species dlffiisivity (assumed constant and 
equal for all species), Yi the mass fraction of 
species i, and v is the gas velocity. These 
equations are to integrated subject to the 
conditions that T = T+ and Yj » Y| + at the 


droplet surface (r = R) as well as T = T„ and 
Yj = Y|„ at infinity. The solutions to these 
equations can be written as 

T-T„ ^ Yj-Y,, 1 - exp(-4»/r) 

T+ - T„ Yu - Y|„ i . exp(-4>/R) 

where the assumption of unity Lewis number 

has been used and $ - r 2 pv(c p /X). By 
applying the conditions for conservation of 
energy and species at J.he droplet surface 

(p + v+L = -X(dT/dr)+ and 8 t = Y| + • 

(D./v)(dYi/dr) + ) it may be shown that <|> - 
R fn(l + B T ) = R ln(l + Bj) where B T = 

c p (T« - T+)/L and Bj « (Y, + - Y^)/(8i - Y l+ ) 

are transfer numbers, 8j is the mass-flux 
fraction of species i in the gas phase, and L = 

Z8jL| is the heat of vaporization of the liquid 
mixture at the droplet surface. It is noted that 
Bj and Bj are required to have the same 
numerical values. 

From conservation of mass principles, it can 

also be shown that k = 8X ln(l+B)/(pLC p ) 
where B *= Bj= Bj. It may be shown that the 
surface mass fraction of a droplet component 

is calculated using Yj+ - (B8 ; + Y^)/(l + B), 

while e is related to B by 

e«=(DLpLC P )/[Un(l+B)]. (6) 


Gas/Liquid Coupling Analysis 

To complete the formulation, the coupling 
between the gas and liquid phases must be 
considered. For simplicity, in this analysis 
phase equilibrium and ideal solution behavior 
will be assumed to describe vaporization. 
Mole fractions (X) and mass fractions are 

related by the identity Xj = OVW i)/(EYj/Wj). 
Gas and liquid mole fractions will be related 
using Raoult’s law and the Clausius* 
Clapeyron relation with constant vaporization 
enthalpy, giving 
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X|. = Xi+exp[Li(l/T + - 1/Tbj)/Ri] . (7) 


I 

V’V = 0 (Continuity) 


In Eq. (7), T b| is the boiling temperature of 
component i at the ambient pressure, and R| 
is the gas constant for species i. The phase 
equilibrium relations used here are not 
essential for the development of this model. 
Nonideal mixtures may be modeled, for 
example, with the use of activity coefficients. 

To solve for the spherically-svmmetrlcal 
- evaporation state, we first specify yio- 
this Is done, values for the surface 
composition in the gas phase are assumed 

(subject to LY i+ *= 1). The surface 

temperature is then iterated until E8i = 1 is 
attained for the gas phase. The gas-liquid 
coupling relations are then used to calculate 
the surface composition on the liquid side of 

the gas-liquid interface. The parameter £ is 
then calculated using Ea. (6), and a value for 
yio based upon the values found above is 
finally calculated from Eq. (3). If this value is 
not sufficiently close to the initially-specified 
value of yiQ, the above procedure is repeated 
with new values for the surface composition 
In the gas phase. Calculations for various 
cases are presented by Shaw [15], where it is 

shown that the 0(e) correction to present in 
Eq.(3) is dominant when y 10 is sufficiently 
small. 

The spherically-sym metrical model has now 
- been completely formulated. The results 
presented above are required for the 
perturbation analysis developed in the next 
section. 


LINEARIZED STABILITY EQUATIONS 
Liquid Phase 

The conservation equations for the liquid are 
assumed to be the incompressible constant- 
property relations appropriate for viscous 
fluids with negligible viscous dissipation. 
These equations are listed below. 


dx ;= - Vp + iiAv (Momentum) 

dl * 

dl-ajAT (Energy) 
dt 


^•cDjAyi (Species) 
dl 


It is to be noted that V represents the gradient 
operator while A represents the Laplacian 
operator. In addition, v is the velocity vector, 

p is pressure, T is temperature, and oti is the 
thermal dlffuslvliy of the liquid. In addition, 

4 = 4 + v-V is the Lagranglan derivative, 
dt dl . , 

Since the droplet surface is regressing, it is 
convenient (as outlined above) to define the 
normalized radial coordinate x = r/R. When 
this is done, the Lagrangian derivative is 
transformed to the form 


<Jt 5t R R dl dx 


where V* is the dimensionless gradient 
operator expressed in terms of x. 

We will now perturb the liquid-phase 
variables by assuming that v = v 0 + v\ 
p-po + p\ t = T 0 + T, yt = yk> + yi’ 
where the subscript "o" denotes the solution 
for the spherlcally-symmetrical state and the 
superscript "prime" denotes the perturbation 

quantity. It Is assumed that v c = 0, p 0 = 
constant, T 0 = constant, and that yi 0 is 
predicted by Eq. (4). The assumption of 
constant T 0 is consistent with the idea that 
early transients associated with initial droplet 
heat-up have decayed. 

Substitution of the perturbed variables 
presented above into the liquid-phase 
governing equations produces evolution 
equations for the perturbation variables (note 
that “primes" will be omitted for brevity). In 
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deriving these equations, the dimensionless 

variables V^ = (8R/k)v', Pl ~ 8p'R 2 /(kHL) 
and Tl » T'/T- are introduced, where T*. is 
the ambient temperature find the subscript "L" 
denotes liquid. The perturbation equations are 
linearized by neglecting products of the 
perturbation quantities or their derivatives, 
and terms involving the spherically* 
symmetrical base solutions cancel. In 

addition, the dimensionless time x defined 
earlier is introduced, though x is rescaled as 

described below (note that exp(-x) = R/Ro). 
The timescale for liquid-phase species 
diffusion is generally much greater than that 
for dlffusioh of energy. A characteristic time 

Xo for species diffusion across the thickness £ 

of the boundary layer is simply x« *= t (on the 

x timescale). For analysis of the equations, 
we therefore introduce the rescaled time 
x=x/e to produce the final form of the 
equations, which are shown below. 

V*>V L = 0 (Continuity) (8) 


A*V l -V*Pl« 0 (Momentum) 

+ « e Le A*T l (Energy) (10) 

e dx dx 


. 1-^. + x^- + "V*V yio 
e dx ™ (ll) 

«eA*yt (Species) 

In deriving Eqs. (8) - (11), it has been 
assumed that the velocity field in the liquid 
phase will be a "creeping” flow (i.e., 
dominated by viscosity) such that fluid 
acceleration can be neglected in the first 
approximation. The variable Le is the liquid- 
phase Lewis number; Le is assumed to be 

large relative to unity and the product £Le is 
treated as 0(1). 

Solutions io Eqs. (8) - (11) will be sought 
assuming that each scalar perturbation 


quantity (Z(x,0,<ti)) can be expressed in terms 
of a normal mode decomposition based upon 
expansions in terms of spherical harmonics 

** n 

Z=£ £ Zn(x)Y7(0,<}>)eP- ? (12) 

n im 0 *n ■ -n 

where Zn depends only on the transformed 

radial coordinate x and Yn*(0,<|>) is the 
spherical harmonic function. The function Zn 
represents a spatial amplitude of the 
perturbation mode h. A normal-mode 
analysis is appropriate here since the 
unperturbed solution is quasisteady in the 

time variable x. A final goal of the analysis 

will be be to calculate the growth factor p n , 
which will determine whether a perturbation 
will grow, decay, or remain steady. 


If we multiply the radial component of the 
momentum equation (V*) by x and then twice 
take the curl, we obtain the equation 



where Vl* is the amplitude of the amplitude 
of the velocity perturbation in the x-direclion 
for the n* mode (for brevity the subscript "n" 
is dropped from all equations involving 
amplitudes). In deriving this equation, use 
was made of the identity 

A[f(x)Y£(0, <{>)]*= 

yni(0 A) + 2.dL . fl 

1 '^Idx 1 Xdx X 2 1 

which is applicable for spherical harmonics 
and where f(x) is an arbitrary function of x. 
Equation (13) is to be solved subject to the 
condition that the solution is finite at x « 0. 
Integration yields 

Vlx * aix nl + aix n+1 (Id) 

where a; and a 2 are integration constants that 
will be determined by considering the 
solutions at the gas-liquid interface (x «= 1). 
Interface conditions will be discussed later. 
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The solution to the energy equation (Eq. 
(10)) may be developed by noting that we 

may reasonably expect that e « 1. To 
maintain a balance between dominant terms In 
Eq. (10), It can be assumed that time 
derivatives are at most 0(1 /Le) In magnitude, 
which implies that the perturbed temperature 
field may be considered quasisteady in the 
first approximation. We will thus drop the 
transient term from Eq. (10), allowing us to 
focus upon the effects of species 
perturbations on stability. Effects of 
temperature perturbations were studied by 
Higuera and Lifldn [16], who found that 
instabilities from temperature perturbations 
were important only when the base solution 
exhibits significant temperature gradients at 
the droplet surface; spatial temperature 
gradients are assumed to have decayed to 
negligible levels for the conditions under 
study here. If we now assume a solution of 
the form of Eq. (12), the resulting ordinary 
differential equation for the amplitude Tlx is 

£Eu + (Z.^]<iXL2s. « 

dx 2 U e Le / dx 

n(n+l)n, n (15) 

■ , • Tlx - 0 . 

X 2 

The point x = 0 of Eq. (15) Is a regular 
singular point, and it is possible to develop a 
Frobenlus series solution about x = 0. When 
this is done, it Is found that the Infinite-series 
solution actually converges very rapidly, and 
as a result the truncated expansion 

Tu = aj[x" + X"* 2 ] (16) 

may be used with little loss in accuracy (as is 
a constant of integration). Applicability of 

Eq. (16) is based upon the condition that eLc 
is not small relative to unity; this condition is 
typically well-satisfied for droplet 
vaporization problems. Though not shown 
here, numerical solutions of Eq. (15) 
compare very well with Eq. (16) for the 
assumed conditions. In developing Eq. (16), 
the condition Tlx *= 0 was applied at the 
droplet center. 


Substituting Eq. (12) into Eq. (1 1) yields the 
evolution equation for the amplitude yiu of 
the liquid-phase species perturbation 


Pn yiU+E X -fa- 
+ (aix n * 1 +a2x n+1 ) (yi.-yio) e (x * 1 ^ 


(17) 


In developing Adequation, the leading-order 
solution for droplet species profiles was used 
(see Eq. (4)). An asymptotic solution to Eq. 

(17) will be sought assuming that e « 1, 
which is generally met for droplet gasification 
in hot environments. Examination of Eq. (19) 
indicates that a boundary-layer structure may 
be expected. The boundary layer will exist 
close to the droplet surface, and it is found 
that an appropriate scaling to allow analysis 
of Ac boundary-layer is provided by the 

transformation x = 1 - ez, where z <; 0(1) In 
Ae boundary layer. Using Ais transformation 
and retaining only Ae largest terms provides 
the leading-order inner problem 


(ai + a 2 ) (yi. - yio) exp(-z) . 


(18) 


The solution to Eq. (18) may be written as 
yiLx « a4texp(s4(l - x)/e) + a 5 iexp($ 5 (l - 

x)/e) - (aj + a 2 )(yi- - yio) exp[(x - l)/e]/p n . 
where 841 and asi are constants and S4 *= [-1 - 
(1 + 4p 0 ) ,/2 ]/2 and s 5 = [-1 + (1 + 

4P„) t/2 ]/2. The leading-order outer solution 
to Eq. (17) is simply yj Lx = 0. To match Ae 
inner and outer solutions we set asi «= 0, 

because for Re(p„) > 0 Ae inner solution 
would grow exponentially fast for increasing 
z. The solution to Eq. (20) may Acrefore be 
expressed as 


yiLx = a 4 iexp[s 4 (l - x)/e] - 

(a! + a 2 ) (yi- - yio) cxp[(x - l)/e]/p n • 


(19) 
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Gas Phase 

For analysis of the gas phase, we will 
introduce the simplifying assumptions that 
density changes, which are a result of 
temperature changes, may be neglected in the 
perturbation equations. We wilfalso assume 

that the vaporization is slow such that <{>/R = 

In (1 + B) « 1 (but where e «= 8 Pi/k is also 
small relative to unity). This allows 
simplified analytical solutions to be found for 
the gas-phase perturbation variables. For 

analysis of the case where ^/R >. 1 and where 
density perturbations cannot be neglected, 
numerical integration of the perturbation 
equations is needed; this analysis is now in 
progress and will be presented elsewhere. 

When $/R « 1, it is easily shown that a 
Reynolds number based upon the gas 
velocity and gas properties at the droplet 
surface Is small relative to unity. Under these 
conditions, we may neglect the convective 
acceleration terms in the momentum equation. 
We will perturb the 6 pherically-symmetrical 

solution by assuming that v = v 0 + v', 

P-Po + P', T = T 0 + r, and Y|*Y lo + Y,' 
where the subscript "o" denotes the solution 
for tire spherically-symmetrlcal state and the 
superscript "prime" denotes the perturbation 
quantity, if we further let x *= r/R and define 
the dimensionless perturbed variables 

V 0 = (po/pU ( 8 R/k)v\ T a = TYT... Y i0 - 

Y|' and Pl = 8 p'R 2 (po/pLy(hM-L)> the 
linearized conservation equations can be 
written as 

V**Vq“ 0 (Continuity) (20) 

A*Vo-V*Po = 0 (Momentum) ( 21 ) 

A*Tq = C(<m) (Energy) (22) 


A Yio *=o(<»/r) (Species) (23) 

where the subscript "G" denotes the gas 
phase. Solutions to Eqs. (20) * (23) can be 
developed in the same fashion as was done 
for the liquid-phase equations. The solutions 
for the amplitudes of the n lh mode of each 
perturbaUon variable are listed below, where 
the condition that perturbations should vanish 

as x — » °o was applied and also where terms 
of order <J/R or smaller were neglected. 

V()x = bi*"" + b 2 x-"- 2 (24) 

Tox = b 3 x-1+O0}>/R) (25) 

Y U x = b 4l x-»-i+0(<|>/R) (26) 

The amplitude in Eq. (24) is for the radial 
velocity component. In addition, b 1 .b 2 .b 3 
and b 4 i are integration constants. 


INTERFACE CONDITIONS 

The solutions for the perturbation quantities 
contain nine unknown constants, which can 
be determined by imposing appropriate jump 
conditions at the liquid-gas interface. Written 
in dimensionless form, equations describing 
the interface conditions in terms of the 
perturbation amplitudes are listed below. 

Vlx — V 0 x (27) 

£C^«J«g&( xi vaj (28) 

-*£- + 2 -d --2 + n(n + 1 ) = 

dx 2 dx J 

Y(£ + 2 £ - 2 + n(n + 1) V u (29) 

+ Mjn(n + 1 )Tl* + My n(n + l)yiLx 

Tlx = To* (30) 
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Vox 


|-(y. 1J ,-e^ i +yi. v u)= 

^-(Yiox + YiooVox)- 


(31) 


(32) 


|-[yau * e + (1 - yi-)Vu] *= 

g*( Y 20x + YmoVqx) - <33) 

yi.L,T_ 

yiu“Pi Y iox*— tt Tl ( 34 ) 

RiT 0+ 2 

(l-y!.)L2T w 

yau - r 2 Y 20 x - ' ?- -- t l (35) 

RaT 0+ 2 

Equation (27) describes overall mass ^ 
conservation across the interface, while Eq. 

(28) equates tangential velocities. Equation 

(29) describes the balance between the gas 
and liquid shear stresses and the surface 
tension gradients that occur from temperature 
and species gradients where Mt = 

8Rai/(T^kllL) and My ■ SRayj/fkpL) arc 

Marangoni numbers and y - (Po/HlXPl/Po). 
Equations (30) equates gas and liquid 
temperatures, and Eq. (31) is an energy 
balance across the interface. Equations (32) 
and (33) are the interface species 
conservation equations, while Eqs. (34) and 
(35) are the linearized versions of the 

equilibrium condition (Eq. (7)) where F j «= 
exp[Li(l/T 0 * - 1 /T m )/Rj]W 0 AV l and T 2 = 
exp[L 2 (l/T 0+ - l/Tb 2 )/R 2 ]WoAV L where Li 
and L 2 , Ri and R 2 , and Tbi and Tb 2 are the 
enthalpies of vaporization, gas constants and 
boiling temperatures of species 1 and 2, 
respectively, and Wq and Wl are average 
molecular weights for the gas and liquid 
phases. In deriving Eqs. (34) and (35), 
perturbations in gas and liquid average 
molecular weights were neglected. 


Substituting the gas and liquid solutions into 
the interface relations produces nine 
equations relating the nine unknown 
constants. The equations can be expressed in 
matrix form as [A][B] — [0], where [B] — 
[at, a 2 » a^ a4 1 • bi, bj, b 3 , b4i, b42] is a 
column vector and [A] is shown on the last 

page of this paper. In this matrix, flj = (n - 
l)/n + My(yj. - yio)/(2p fl ). = (n + 2)/(n + 

1) + My(yj. - yio)/(2Pn), fi 3 = l + 

[16eLe(2n + 3)]->, - (KJXq) ( n“? (n + 

2) /[16eLe(2n + 3))}, \j/i » yi. LiT«/(RiT 0t 2 ) 
and V 2 “ (1 - yi-)L< 2 T«/(R 2 To+ 2 ). 


STABILITY PREDICTIONS/DISCUSSION 


The stability of two-component droplet 
gasification can be evaluated by setting the 
determinant of [A] equal to zero and solving 


for the growth factor p n . Because of the large 
number of terms in the matrix, this was 
accomplished using a symbolic calculation 

software package 122]. The variable P„ was 


found to depend upon a large number of 
terms (roughly 100), and for brevity the 
dispersion relation is not listed here. 


Surprisingly, it was found that the dispersion 
relation did not depend upon either of the 
Marangoni numbers defined earlier (M A and 
My); they cancelled out of all calculations in 
the determinant. This was verified directly by 
calculating the minors associated with the 

terms flj, flj, -M A (V2 and -My/2 in [A]; 
the minors for -M A H3/2 and -My/2 were 
Identically zero, while the minors for and 

were exactly equal. In addition, the 
dispersion relation predicted that droplets 

should be absolutely stable, that is, p n < 0 in 
general. In retrospect, this result might have 
been expected. To clarify why this is so, we 
should recall that the leading-order solution 
(yi = yio + (yi. - yio)oxp(-z)) was used for the 
unperturbed state, while Eq. (3) shows that 
the leading-order expression for the mass- 
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flux fraction of species i leaving the droplet is 

simply 8i = yio, i.e., 8i is dependent only 
upon the initial mass fraction of species i and 
is independent of variations in other 
variables. As a result, in this diffusion- 
controlled limit for liquid species transport, 
small perturbations in liquid-phase velocity 
and species profiles do not influence the 
mass-flux fraction of a species leaving a 

droplet, and 6i will be unchanged over the 
surface of a droplet as will the liquid-phase 
composition at the surface. By Eq. (5), the 
gas-phase solution will therefore be 
unchanged, and small perturbations in the 
liquid phase will not significantly alter gas- 
phase species or temperature fields so long as 

8j « yio is an adequate approximation. Under 
these conditions, all perturbations die out 
from viscous decay of velocity fields. 

It should be noted, however, that these 
conclusions are not expected to be valid for 
all cases. For example, it is possible for the 

0(e) correction in Eq. (3) to be important, 
and when this is so the analysis presented 
here may not be valid. This can happen, for 
example, when species i is initially present in 
very small amounts and is of low volatility 
relative to the other liquid species so that it 
must build to appreciable levels at the droplet 
surface before it can vaporize. Another 
possibility is when species i was not even 
present initially in the droplet, but was first 
absorbed and then gasified from the droplet 
surface (e.g., water condensation onto 
alcohol droplets). For analysis of the stability 
characteristics of droplets under these 
conditions, further analysis is necessary. 
These analyses are being performed, and will 
be published at a later date. 

Even though these results were derived for 

slow gasification rates (<}>/R « 1), it is 
expected that they will be valid for rapid 

gasification (<Ji/R £ 1) as well as for droplet 
combustion. So long as the liquid species 

transport is diffusion controlled and 8] - yio 
is a good approximation for the liquid phase, 
small perturbations in liquid-phase variables 
(e.g., velocities) will not significantly affect 


8|. As a result, small perturbations in liquid- 
phase . variables will not induce surface 
temperature or species gradients, and velocity 
perturbations will decay from viscous effects. 
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